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ABSTRACT 

We compare cluster scaling relations published for three different samples selected via X-ray and 
Sunyaev-Zel'dovich (SZ) signatures. We find tensions driven mainly by two factors: i) systematic 
differences in the X-ray cluster observables used to derive the scaling relations, and ii) uncertainty in 
the modeling of how the gas mass of galaxy clusters scales with total mass. All scaling relations are 
in agreement after accounting for these two effects. We describe a multivariate scaling model that 
enables a fully self-consistent treatment of multiple observational catalogs in the presence of property 
covariancc, and apply this formalism when interpreting published results. The corrections due to 
scatter and observable covariance can be significant. For instance, our predicted Ygz~Lx scaling 
relation differs from that derived using the naive "plug in" method by « 25%. Finally, we test the 
mass normalization for each of the X-ray data sets we consider by applying a space density consistency 
test: we compare the observed REFLEX luminosity function to expectations from published Lx~M 
relations convolved with the mass function for a WMAP7 flat ACDM model. 
Subject headings: cosmology: clusters 



1. INTRODUCTION 

The study of scaling relations of galaxy clusters is 
one of significant importance, both from a cosmolog- 
ical and an astrophysical point of view. Cosmolog- 
ically, mass-observable scaling relations are a funda- 
mental component of all work that exploits the abun- 
dance of galaxy cl usters for constraining cosmologi- 
cal parameters (e.g. iHenrv et al.l 120091: [Vikhlinin et al.l 
l2009atlMantz et al.li2010aHRozo et al. 1120101 ) Astrophys- 
ically, the existence o f a self-similar model (|Kaiserlll986t 
IBohringer et al.ll2012[ ) allows one to use departures from 
this self-similar expectation as a probe of feedback and 
non- thermal processes in cluster and galaxy formation 
(e.g.lRowlev et al.ll2qol iMagliocchetti fe Bruggenll2007l: 
Mantz et al l l2010bl: iMittal et al.l 120111: iMaughan et al l 
2012t lEckmiller et all 120111 ). Likewise, the scatter of 
galaxy clusters about the scaling relations can also be 
a direct probe of the state of the intra-cluster gas (e.g. 
the presence or absence of a cool c ore) and/or dynam- 
ical state (e.g. iFabian et al.l H99l lHartlev et al.l 120081: 
iRasia et al.1 120111: iKrause et al.l 120 lW . Indeed, these 
ideas have spurred large suites of simulations in which 
the physics contributing to the state of the intra-cluster 
gas are systematically varied, both to give a range of 
plausible evolution models, and to help guide theoret- 
ical interpretations of current observations (e.g. iNagail 
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All this work, however, is critically dependent on our 
ability to robustly measure empirical cluster scaling re- 
lations. Here, we take a pragmatic approach to esti- 
mating the level of systematic differences in observed X- 
ray and SZ scaling relation by comparing the published 
values of three differe nt sets of works: a Planck and 
XMM based analysis dPratt et all 120091 : lArnaud et all 
120101 : iPlanck Collaboration! 120 llbl ). henceforth referred 
sim ply as Pll-LS, a Chandra and Planck b ased analy- 
sis (jVikhlinin et al.ll2009aHRozo et al.ll2012al) . henceforth 
referred to as V09, and a second, i ndependent Chandra 
based analysis (jMantz et al.l 120 lObl) . henceforth referred 
to as M10. As is demonstrated below, these three works 
are often in tension with one another. Understanding the 
origins of these tensions is of paramount importance for 
validating the use of galaxy clusters as probes of precision 
cosmology. 

This is the second in a series of papers that aims to 
perform a detailed study of local cluster scaling rela- 
tions from X-ray a nd SZ selec t ed clus ters catalogs. In 
the first paper — iRozo et"aLl (|2012bD . henceforth pa- 
per I — we demonstrate significant systematic differences 
between different groups on raw cluster observables, in- 
cluding X-ray luminosity, Lx, temperature, Tx, and gas 
thermal energy, Yx- More importantly, hydrostatic mass 
estimates vary substantially. 

Paper I characterizes relative offsets in these quanti- 
ties, and that information forms the foundation for the 
current work. Indeed, we demonstrate below that the 
tension in the scaling relations between the three works 
we consider is ultimately sourced by the systematic dif- 
ferences in cluster observables identified in paper I. 

Having fully understood the origin of the tension be- 
tween the three groups, we then proceed to derive a 
self-consistent set of multi-variat e scaling relations from 
each of these data sets (see also lWhite et al.|[2010|) . By 
"self-consistent", we mean that our propagated seal- 
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ing relations are explicitly derived from a probabilistic 
model that account for the scatter and possible covari- 
ance between clusters obscrvablcs, a rising for instance 
from local large scale structure (e.g. IWhite et al.l 120101 : 
iNoh fc Cohr][20TTl [20ll lAngulo et al.ll2012j ) Most work 
to date relies on simple "plug-in" methods to propa- 
gate scaling relations. That is, given two scaling rela- 
tions Y(X) and Z(Y), the Z-X scaling relation is as- 
sumed to take the form Z(Y(X)). Recent work with nu- 
merical simulations and/or Monte Carlo analyses have 
made it clear, however, that this naive method is gener- 
ically biased, with the error depending on both the 
scatter and covar i ance between the observables at hand 
(IRozo et alj|2009t IWhite et all 1201(1 iBiesiadzinski et al.l 
l2012tlAngulo et al.ll2012D . The two specific cases we con- 
sider are using the M—Yx and Ygz-Yx scaling relations 
to derive the Y$z~M scaling relation, and using the Lx~ 
M and Y$z~M scaling relations to derive the Ygz~Lx 
scaling relation. The latter example is one where ignor- 
ing scatter corrections can result in biases as large as 
20% - 30%. 

Having derived the cluster scaling relations for each 
of the three data sets we consider, and having identified 
the origin of the differences them, we turn to investigate 
which of these analyses, if any, is consistent with cos- 
mological expectations, for a WMAP7 +BAO+ff flat 
ACDM best fit cosmological model of iKomatsu et al.l 
(|2011l ). Specifically, we consider each of the published 
Ly-M scaling relat ion in turn, convolve them with the 
iTinker et all (|2008l ) mass function appropriate for the 
afformcntioned cosmology, and then compare the result- 
ing predicted abundance with the X-ray luminosity func- 
tion from the REXCESS catalog (jBohringer et all 120021 
l200l . 

In addition to being interested in the X-ray and 
SZ scaling relations of galaxy clusters in their own 
right, one of our main motivations for pursuing this 
work is to take a closer loo k at the recent results by 
iPlanck Collaboration! (|2011cD . In that work, it was found 
that the Ysz-A^nn scaling relation of o ptically-identified 
maxBCG clusters (jKoester et alJl20"07t ) were inconsistent 
with predictions d erived from a combina tion of optical 
(IRozo et al.l l2009f) and X-ray data sets ([Arnaud et al.l 
120101) . Because these predictions depend on the X-ray 
and SZ scaling relations of galaxy clusters, it is impera- 
tive that systematic uncertainties associated with these 
predictions be quantified, and that the scaling relations 
be self-consistently propagated in the presence of non- 
zero scatter. The results from this comparison will be 
presented in a subsequent paper, Rozo et al. (in prepa- 
ration), henceforth referred to simply as paper III. For 
the purposes of this work, this means we will focus our 
attention on cluster scaling relations at z = 0.23 exclu- 
sively, the median redshift of maxBCG clusters. 

The layout of the paper is as follows. In § [5] we de- 
scribe the data used in this work. This data takes the 
form of the Lx~M, M-Yx, and Ysz-^x scaling relations 
derived from each of the data sets we consider. Section 
[3] compares these three scaling relations to each other, 
and demonstrates that after correcting for the system- 
atic offsets in X-ray obscrvablcs identified in paper I — 
and tilting the M10 relations to an / gas oc M 015 model 
— all scaling relations are in good agreement with each 
other. We then turn to self-consistently deriving the 



Ysz~M (section @| and Ygz~Lx (section [SJ scaling re- 
lations from the input scaling relations summarized in 
section [2] Section [5] also explicitly compares our self- 
consistently derived s c aling r elation with the data from 
IPlanck Collaboration! (|2011bf ). Section |6] tests whether 
the scaling relations from each of the three works are con- 
sistent with the observed cluster luminosity function in 
the current best-fit ACDM cosmology. Section[7]presents 
a summary of our results. 

In all cases, scaling relation parameters are computed 
assuming a flat ACDM cosmology with D, m = 0.3 and 
h = 0.7. When comparing abundances to cosmolog- 
ical expectations for a WMAP7 +BAO+j?o cosmology , 
we use the best fit model from IKomatsu et al.l ([201 ID 
- which has Q m = 0.275, h = 0.702, n s = 0.968, 
and as = 0.816 — to compute the predicted halo mass 
function. As discussed in section [6j the scaling rela- 
tions we use are still those that were computed assuming 
O m = 0.3. This is slightly inconsistent, but is necessary 
in the absence of precise knowledge of the degeneracies 
between each set of scaling relation parameters and cos- 
mological parameters. We note, however, that because 
the changes in cosmological distances to low redshifts are 
very mild between the two cosmological models, we ex- 
pect this inconsistency will only impact our predictions 
at the few percent level at most. This level of uncertainty 
is not enough to modify our conclusions. 

Unless otherwise noted, all total masses, M, employ 
the M500 convention, where M500 is defined as the mass 
within a radius, -R500, that encompasses a mean interior 
density of 500 times the critical density of the universe, 
p c (z) = 3H 2 (z)/8ttG. 

2. DATA 

We consider the cluster scaling relations as measured 
by three groups. The first is the Pll-LS data set, which 
employs Planck and XMM data, and is comprised of 
IPlanck Collaboration! (|2011b[ ). and th e associated pub- 
lications based on R EXCESS clusters ()Pratt et al.ll2009t 
lArnaud et al.ll2010( ) . We note that in paper I we demon- 
strated tjiat : the_£_<_0JJ_ajid z g [0. 13, 0.3] galaxy clus- 
ters in IPlanck Collaboration! (|2011bf ) appear to be sys- 
tematically difference. Consequently, and motivated by 
our ultimate goal of investigating the maxBCG Ysz-A^oo 
scaling relation, we will also be be explicitly consider- 
ing the subset of Pll-LS galaxy clusters in the redshift 
range z € [0.13,0.3], as appropriate for maxBCG sys- 
tems. We refer to this subset of the Pll-LS data set as 
P ll-LS(z=0.23). T he second d ata set, V09, is comp rised 
of lVikhlinin et all (j2009al) and IRozo et~aTI (|2012af ). but 
relies on Chandra rather than XM M data. Finally, w e 
consider the X-ray analysis of M10 (jMantz et al.ll2010bT ). 
which is also Chandra based. Because this is the sec- 
ond in a series of papers, we will simply refer the reader 
to paper I for a more detailed description of each of the 
cluster data sets on which the scaling relations analysis is 
based. Here, we focus exclusively on the scaling relations 
reported in these works. 

Table [1] summarizes the Lx~M, M-Yx, and Ysz~Yx 
scaling relations as quoted in each of the above works, 
though modified when necessary to match the definitions 
for cluster observables adopted in this work (see below). 
We always adopt the pivot points reported by each indi- 
vidual work. The only exception to this rule is that in 
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TABLE 1 

Input Cluster Scaling Relations at z = 0.23 



Relation 


X0 


a 


a 


°lni/i|x 


Citation 


Data Set 




4.8 


1 1 r 4- n nq 


1.61 ± 0.14 


n Qofi + n 039 


Vikhlinin pt al f9n09a > l 

V Ixvlllllllll t_ L ell. \ — V )\j iJ<\ J 


V09 


Z/X-.M500 


2.0 


0.08 ± 0.08 


1.62 ±0.11 


0.411 ± 0.070 


Pratt et al. (2009^" 


Pll-LS 


Lx-Af500 


10.0 


2.11 ±0.18 


1.34 ±0.05 


0.414 ± 0.044 


Mantz et al. f^oTob) 


M10 


Lx-M 500 


4.0 


0.98 


1.52 




Reference 




M 5Q0 -Y X 


3.0 


1.53 ±0.04 


0.57 ±0.03 


< 0.07 b 


Viklilinin et al. (2009a) 


V09 


M 50 o-Y x 


2.0 


1.23 ±0.02 


0.56 ±0.02 


< 0.09~ 


Arnaud et al. (2010) 


Pll-LS 


M 500 -Y x 


10.0 


2.25 ± 0.12 


0.68 ± 0.04 


0.072 ±0.011 


Mantz et al. (2010b) 


M10 




4.0 


1.65 


0.6 




Reference 




D'iYsz-CYx 


8.0 


1.877 ±0.028 


0.916 ± 0.035 


0.082 ±0.035 


Rozo ct al. (2012a) 


V09 


D$Y SZ -CY X 


10.0 


2.341 ±0.038 


0.828 ± 0.057 


0.167 ±0.039 


Rozo et al. (2012a) 


Pll-LS(z=0.23) 


B%Y SZ -CY X 


10.0 


2.100 ±0.09 


1.0 


< 0.15^ 


This work 


M10 


d|y sz -CY x 


10.0 


2.303 


1.0 




Reference 





a In all cases, we assume the ij)— x relation takes the form {\ntp) = a + a ln(x/xo)- Our choice of units are 10 14 Mq 
for mass, 10 44 ergs/s for L x , 10 14 M0keV for Yx, and 10 -5 Mpc 2 for D^Ygz an d CY X . Unless otherwise noted, 
we set xo to the reference scale in the cited work. All scaling relations are evaluated at z = 0.23, the median 
redshift of the maxBCG cluster sample. 

^ 1 Vikhl inin et al. (2009a|) only state that the scatter is unde tectable given the error s on hydrostatic mass estimates, 
but that this is consistent with 7% scatter as predicted by Kravtsov et al. (2006). We implement this scatter in 
o ur analysis as a unifo rm prior in the variance with the maximum value quoted above. 

c Arnaud ct al. (2007) quote a scatter of 0.087, but provide no error bars. We implement this scatter in our 
analysis as a uniform prior on the variance using the maximum value quoted above. 

d Uncertainty in the scatter is implemented as a uniform prior on the var iance with the maximum valu e quoted 
above. The maxim um value is chosen to be close to that derived from the Planck Collaboration (2011b) data by 
IRozo et al] H2012af) . 

the Lx~M relation from V09, whe r e we s et Mq to the 
median mass of the lVikhlinin et all (|2009al) low redshift 
cluster sample, since the published relation used 1 M Q 
as its pivot point. We emphasize that in all cases, given 
two observables Y and X, our definition of the Y—X scal- 
ing relation is the expectation value (\nY\X). So, for 
instance, in the case of the Lx~M relation derived from 
X-ray samples, we only ever consider the Malmquist-bias 
corrected relations. This uniformity of definition is cru- 
cial for a self-consistent analysis of multi-variate cluster 
scaling relations. When considering stacked relations, 
we also take care to correct for the expected difference 
between hi(Y\X) and (lnY|X) within the context of a 
log-normal scatter model. 

Because our final goal (paper III) is to use the 
results from this work to investigate the discrep- 
ancy between theory and o bservations uncovered by 
IPlanck Collaboration! ()2011bl ). we have evaluated all 
scaling relations at z = 0.23, the median redshift of 
the maxBCG clus ter sample. The Y^z-Yx relation in 
IRozo etldl (|2012al ) is constrained only at z w 0.1, and 
we assume no redshift evolution when extending this 
result to z = 0.23, consistent with self-similar evolu- 
tion. We note that all three works have different evo- 
lution terms in their scaling relations. For instance, 
for £ y-M. iPratt et al l (120091) as sume L x oc E(z) 7 ^ 3 
while I Vikhlinin et all (|2009aD and iMantz et all (|2010bf ) 
find a best fit evolution Lx oc E{z) 1 - 85±0A2 and Lx oc 
-E(z) 2 ' 34±a05 respectively. The self-similar expectation 
for soft-band X-ray luminosities is approximately Lx oc 
E(z) 2 . The evolution factors of Pll-LS and M10 are 
nearly identical, with the evolution relative to V09 scal- 
ing as E(z) 0A9 ~ 1.06. That is, the relative evolution 
from z = to z = 0.23 can induce up to a 6% differ- 
ence between the various works. This is also an extreme 
case, since all scaling relations have a redshift pivot point 
z > 0. Because these ~ 5% differences are typically much 



smaller than that the differences between works, we will 
simply ignore them here, evaluating all scaling relations 
at z = 0.23 using each group's E(z) evolution factors. 
Determining what the correct evolution factor is for each 
scaling relation is beyond the scope of this work. 

We have also homogenized all data to common defini- 
tions: Lx is defined as the total X-ray luminosity within 
i?50o in the [0.1, 2.4] keV band, Tx is defined as the spec- 
troscopic temperature in an an annulus R G [0.15, l]i?5oo, 
M gas is the total gas mass within i?soo, Yx = M gas x Tx, 
and Ysz is the integrated Compton parameter within a 
sphere of radius i?soo- Whenever our definitions do not 
match the definitions employed in any of the works we 
consider, all data is rescaled as described in paper I. The 
data in Table [1] — in conjunction with the systematic 
offsets we estimated in paper I — form the basis of our 
analysis. 

The different groups rely on different analysis tech- 
niques and different mass proxies when estimating cluster 
scaling relations. For the Lx~M analysis, V09 and Pll- 
LS (Via IPratt etaLl [2009] ) rely on Y x as a mass proxy, 
while Mantz relies on M gas . In all cases, the mass-proxy 
(Yx or M gas ) relation with mass is calibrated using hy- 
drostatic mass estimates of relaxed galaxy clusters. V09 
employs a likelihood fitting method that explicitly incor- 
porates Malmquist bias, and tests their fitting routine 
using Monte Carlo cluster samples. The M10 analysis is 
very similar in spirit to the V09 analysis, but also ex- 
plicitly incorporates the cosmological dependence of the 
cluster abund ance function in t he cluster scaling rela- 
tions. Finally. |Pratt et aLl (|2009D quote two different fits, 
corresponding to BCES(Y\X) and SCE'S'-orthogonal 
fits. BCES(Y\X) is closest to a likelihood fit in the 
absence of Malmquist bias. Since this fit is the one that 
is most comparable to those of V09 and M10, we only 
employ the BCES (Y\X) Malmquist-bias corrected fit of 
IPratt et al.l ()2009l ) in this analysis. We caution that this 
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TABLE 2 

Mean Log Differences in X-ray Properties for Sample Pairs 



Property 


M10-V09 


P11-LS-V09 


P11-LS-M10 
Low z 


P11-LS-M10 
High z 


P05-V09 


Tx 

Yx ab 


0.12 ±0.02 
0.03 ±0.02 

-0.03 ± 0.02 
0.11 ±0.11 


-0.01 ±0.02 
0.03 ±0.02 
-0.13 ±0.02 
-0.15 ±0.03 
-0.12 ±0.02 
-0.14 ±0.03 


-0.12 ± 0.02 
-0.02 ± 0.03 

-0.05 ± 0.07 
-0.14 ±0.03 


-0.10 ±0.03 
-0.04 ±0.02 
-0.14 ±0.05 
-0.19 ±0.05 
-0.37 ±0.07 
-0.34 ±0.06 


-0.18 ± 0.05 


a Offset computed after outlier removal. 



b Offsct computed after correction to a common aperture. 
c Relaxed/cool core only. 

d Relative to paper I. the offsets here are corrected by an 11% shift due to an update in 
Chandra calibration. This update was not applied to the masses quoted in the Tables 
in M10, but is applied for the scaling relations. 

c Non-relaxed/no cool core only. 



is still not directly comparable to likelihood based ap- 
proaches, since neglecting intrinsic scatter in the cost 
function being minimized gives cluste rs with small er- 
rors undue weight in the fit (see also lAndreon fc Hurnl 
|20f0j)P| Importantly, as we show below, these differences 
in fitting methodology are not the primary drivers of the 
tension between the various scalin g relation analyses. 

Turning to the M-Yx relation. lArnaud et alt (|2010Q 
rely on hydrostatic mass estimates of relaxed galaxy clus - 
ters to calibrate M-Yx- Because lArnaud et al.l (|2010l ) 
did not quote the corresponding scatter, we rely instead 
on the v alue quoted in th e ir prev ious work, lArnaud et al.1 
(|2007fl . IVikhlinin et~aTl (|2009af ) followed a similar ap- 
proach to determine M-Yx, but also fail to quote a con- 
straint on the scatter, noting only that the data is consis- 
tent with no scatter, a s well as with the predictions from 
IKravtsov et al.l (|2006|) . Thus, we place a n upper limit on 
the sc at ter based on the valu e reported in lKravtsov et al.l 
(|2006Q . iMantz et all (|2010bh uses M gas as a mass proxy, 
and the M-M gas relation is itself calibrated based on hy- 
drostatic mas s estimates o f relaxe d galaxy clusters. We 
note too that IMantz et all (|2010b[ ) report Y x -M rather 
than M-Yx- To obtain M-Yx, we simply invert their 
relation. The correctio ns d ue to scatter are 2% in the 
amplitude (see equation I A5j) . and are completely negligi- 
ble both relative to their quoted uncertainty, and for the 
purposes of this work. 

Finally, turning to the lsz~^x scaling relation, we 
write this as 

(\n{D\Y S7 ) \CYx) = a sz + a sz \xiCYx. (1) 

where Da is the angular diameter distance, and C is a 
normalization constant, 



C 



1 



1.407 x 



10" b Mpc" 



(2) 



m e c" Lim p 10 14 M Q keV 

where ot is the Thompson scattering cross section, m 



for instance, when mea- 
£ [0.12,0.3] clusters in 



8 We see very clear evidence of this, 
suring the Ysz~^x relation from the 

IPlanck Collaboration! COllbl) : a BCES(Y\X) fit results in steeper 
scaling relations than a likelihood fit because a few very massive 
clusters with extremely tight error bars dominate the penalty func- 
tion being minimized. By including intrinsic scatter as a free pa- 
rameter in the fit — as opposed to computing it a posteriori based 
on a fit including only statistical errors — these clusters are de- 
weighted, which flattens the resulting relations. 



is the mass of the election, m p is the mass of a proton, 
and ll = 1.156 is the mean molecular weight of the intra- 
cluster medium for solar abundance (variation in metal- 
licities indues a 1% uncertainty, completely negligible 
for our purposes). 

Our constraints for the Ysz~Yx sc aling relation come 
from two sources. IRozo et "all (|2012aD constrain the Ysz - 
Yx ratio using Planck and Chandra and this is the re- 
sult we use later t o estimate the Y $ z~M r elation in con- 
junction with the IVikhlinin et al.l (]2009aD data. Turn- 
ing to the Pll-LS data set, we noted in paper I that 
the Pll-LS Ysz/Yx ratio is different for z < 0.13 and 
z > 0.13 clusters, s o rathe r than using the fit from 
IPlanck Collaborationl (|2011bl) for the lsz _ ^x relation, 
we have refit the Pll-LS data for galaxy clusters in 
the z g [0 . 13,0.3] redshift range using the method of 
IRozo et all ()2012al ). Our best fit scaling relation for the 
Pll-LS(z=0.23) data is included in Table [TJ 

It is important to emphasize that by choosing to use 
the Ysz-Yx scaling relation derived from the Pll-LS data 
in our predictions for the Y$z~M (and eventually A^oo - 
M ) scaling relations, we are diver ging from the treatment 
of IPlanck Collaborationl (|2011d ) for deriving the Ysz~ 
A^oo scaling rel ation. In that work, the Y$z~M relation 
used was that of lArnaud et al.l (|2010D , which was derived 
from X-ray clusters alone. A more detailed discussion 
about the X-ray exp ectations for the Y sz~M scaling re- 
lation is presented in lRozo et al.l (|2012aD . Here, the point 
that we want to emphasize is that by working with the 
observed Ygz~Yx scaling relation, we are explicitly test- 
ing whether all scaling relations are self-consistent, in the 
sense that a multivariate log-normal model can explain 
all observational data. Discussion of the consistency of 
the Ygz~Yx scaling rela t ion wit h X-ray expectations can 
be found in IRozo et all ()2012al ). 

For completeness, we have also estimated the 
D\Y$z/CYx ra tio that would be der i ved fro m a com- 
bin ation of the Planck Collaborationl (|2011bO data and 
the IMantz et al.1 (|2010b( ) data ( i.e. we are explic tly set- 
ting asz = 1 for the M10 data set). For all galaxy 
clusters shared by the M10 and Pll-LS cluster samples, 
using the formalism from paper I we rescale the M gas 
values in M10 to the -R500 aperture adopted by Pll- 
LS. We then directly compute Yx for these systems, and 
use bootstrap resampling to compute the mean value of 
\n(D A Ysz/CYx) at this aperture, along with the associ- 
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ated uncertainty. We find (\n(Y sz /Y x )) = -0.16 ±0.09. 
We have confirmed that a pplying this exact method to 
the lVikhlinin et all (|2009al ) data results in ln(F sz /Y x ) = 
—0.20 ± 0.04, in perfect agreement with the more de- 
tailed analysis of IRozo et al.l (|2012aD . We do not pro- 
vide a scatter estimate since the errors in IMantz et al.l 
(|2010b|) represent the expected uncertainty in lx, rather 
than just measurement errors. E.g. the error bar in- 
cludes the expected scatter in gas clumping from cluster 
to cluster, which alone contributes 5% scatter, already 
larger than the typical measurement error. The upper 
limit to the scatter quoted in Tab le [J is assumed based 
on the scatter d erived from the iPlanck Collaboration! 
(|2011bf) data by IRozo et al.l (|2012al) . Note, however, 
that comparison to the Pll-LS and V09 values sug- 
gests the quoted upper limit is reasonable. The best 
fit D\Ysz/CYx is 0.82, corresponding to an amplitude 
parameter a = 2.10 for the Ysz-^x relation. 

For future reference, Table[3]shows the offsets in X-ray 
derived cluster observables from paper I. In fact, table 
[2] is a reproduction of Table 3 from paper I, albeit with 
some sli ght modifications. First, because the Lx~M re- 
la tion oflPratt et al.l (l2009j) re lies on the mass calibration 
of IPointecouteau et al.l (|2005h — henceforth P05 — we 
have also explicitly added the mass offset between P05 
and V09 for the purposes of our comparison (the off- 
set was computed in paper I, but not included in Table 
3 of that work). In addition, the scaling relations in 
IMantz et al.l (|2010bf ) were corrected for a Chandra cal- 
ibration update, which uniformly shifts the masses of 
galaxy clusters by ~ 11%. Because the cluster data ta- 
bles in IMantz et ail (|2010bO were not corrected for this 
calibration update, but the scaling relations were, we 
have corrected the observed P11-LS-M10 and M10-V09 
offsets fro m paper I values to r eflect the calibration em- 
ployed by IMantz et al.l (|2010bf) in their scaling relation 
analysis. This calibration update significantly improves 
the agreement between the Pll-LS and M10 masses. Be- 
cause Mantz quotes iJOSVlT-calibratcd luminosities and 
M gas , there is no corresponding shift in Lx or M gas . 
While M gas (i?5oo) does shift due to aperture effects, Ta- 
ble [5] quotes the M gas offset at a fixed aperture, so an 
aperture correction is not necessary. 

3. COMPARISON OF CLUSTER SCALING RELATIONS 

Below, given a scaling relation Y-X, we will compare 
the 68% confidence regions in the Y-X plane correspond- 
ing to the probability distribution P(\nY\\nX). These 
regions are estimated by performing 10 5 Monte Carlo 
draws of the amplitude and slope of each published scal- 
ing relation. We then evaluate (In Y| lnX) along a log- 
arithmically spaced grid in X, and compute the stan- 
dard deviation of the resulting (In Y| In A) values at each 
point. This defines our estimate of the 68% confidence 
region of each cluster scaling relation in the Y-X plane. 

3.1. The Lx -M Relation 

The top panel in Figure Q] compares the Lx-M^qq rela- 
tions derived from the Pll-LS, V09, and M10 data sets. 
We remind the reader that the Lx~M rel ation we as- 
sociat e with the Pll-LS data set is that of lPratt et al.l 
(|2009Il The width of these bands reflects the quoted er- 
rors in each of the works being compared, as summarized 
in Table [U 



To more easily compare these relations, we compute 
the offset of each relative to a mean relation. The ampli- 
tude of the mean relation is defined as the average value 
of InLx, evaluated at mass A/500 = 4 x 10 14 M Q , as esti- 
mated from the Pll-LS, V09, and M10 best fit scalings. 
The slope of the mean relation is likewise defined as the 
mean of the slopes of the three scaling relations. This 
mean scaling relation is only intended to be a reference 
to be subtracted from the three scaling relations: we do 
not ascribe any physical interpretation to it. More specif- 
ically, we do not advocate interpreting the mean scaling 
as more likely to be correct than any of its constituent 
elements. 

The M10 errors of Tablc[T]arc larger in amplitude than 
those quoted for the Pll-LS and V09 data sets because 
they include systematic uncertainties, whereas the V09 
and Pll-LS errors are measurement errors only. The 
published M10 slope, however, has an error of only 0.05, 
compared with 0.11 and 0.14 for Pll-LS and V09, re- 
spectively. These differences in the relative errors of the 
amplitude and slope explains why the Pll-LS and V09 
contours have an hourglass shape, while the M10 con- 
tours are nearly parallel lines. 

The lower panel in Figure Q] shows the 68% confidence 
intervals of the Lx~M relations after subtracting the 
mean relation. We see that the Pll-LS scaling relation is 
parallel to but clearly offset from the V09 scaling relation, 
while the M10 relation is flatter. We can test whether 
the difference in amplitude between these various scaling 
relations are due to differences in the raw X-ray clus- 
ter observables simply by shifting the amplitude of all 
the relations by the systematic offsets quoted in Table 
[2] If A InLx = (Inix)B — ( m ^x)A i s the mean offset 
between samples B and A in log-luminosity, and with 
AlnM defined similarly for total mass, then the scaling 
relation from sample A can be observable-offset corrected 
to sample B by adding the constant [A In Lx — «aA In M] 
to (InLx) a- Assuming that the only difference between 
the two data sets A and B are the observable offsets char- 
acterized in paper I, applying these corrections to every 
data set should bring all scaling relations into good agree- 
ment with each other. For specificity, we observable- 
offset correct every scaling relation to the Pll-LS scaling 
relation. 

These amplitude shifts do not affect the slopes of the 
relations, and therefore cannot explain the tension be- 
tween the slope derived by M10 and the remaining two 
data sets. Instead, the difference in slope can be traced 
to differences in ,f gas - Af between the y a rious works. 
Based o n the results by lAllen et~aT] (|2008D . IMantz et al.l 
(|2010bft assume that / gas = /o, a constant value. In pa- 
per I, we saw that given a different, power-law model, 
/ gas (M) ~ M 7 , the relation between the mass M of that 
model and that estimated by M10 is 

where M p is the pivot mass of the f gas -M relation, and 
7 is the index of the f gas -M power-law. If tp is a cluster 
observable where the M10 scaling relation is tp oc M aM1 ° , 
then the scaling relation the M10 scaling relation for an 
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Fig. 1. — Top panel: 68% confidence band for the L-&—M relation 
for the Pll-LS, V09, and M10, data sets, as labelled. All scaling 
relations are evaluated at the median redshift of the maxBCG sam- 
ple, 2 = 0.23. Bottom panel: As per the top-panel, after subtract- 
ing out a reference mean scaling relation (see text for the definition 
of the mean scaling relation). Tensions are evident among the scal- 
ing relations. 



/g as oc M 1 would be 
tp = A 



M \ 1 ' 677 M 



M r , 



(4) 



Here, Mq, A, and «mio are the pivot poin t, amplitude, 
and slo pe of the tp—M relation reported bv lMantz et al.l 
(|2010bD . The slope a relative to the mass M in this new 
/g as model is 



aAfio(l + 1-677). 



(5) 



In the work of V09 and Pll-LS, the slope 7 of the / ga s _ 
M relation is 7 « 0.1 — 0.2. Setting 7 = 0.15, and given 
otMXQ = 1-34 ± 0.05, our model predicts an Lx~M slope 
a = 1.68 ± 0.05, resulting in an excellent match to the 
values of Pll-LS and V09. Evidently, the reason that 
the M10 scaling relation appears flatter is that the / gas 
model employed by M10 is not consistent with that of 
Pll-LS or V09. 
Figure [2] shows the difference in the Lx~M relation 



Fig. 2. — 68% confidence band for the Ly^—M relation derived 
from the Pll-LS, M10, and V09 data sets, as labelled. We have 
subtracted out a reference scaling relation as per the bottom panel 
of Figure [T] All scaling relations have been displaced along the 
Tx and M axis by the observable offsets relative to the Pll-LS 
and P05 data sets characterized in Paper I. In addition, the M10 
scaling relation has been tilted from a constant / gaa model to a 



model with / gi 



M 0,15 . These corrections improve agreement 



among the scaling relations. 



of V09, M10, and Pll-LS, after accounting for the ob- 
servable offsets in Table [2] The M10 relation is also 
tilted to an / gas oc Af 015 model as described above0 
For M10, we show the Lx~M relation after shifting the 
M10 values using the P11-LS-M10 offset (dashed lines), 
and then also applying the offset be tween Pll-LS and 
P05, which is necessary because the iPratt et al.l (|2009T ) 
scaling relation is calibrated based on P05. With these 
corrections, all scaling relations are in excellent agree- 
ment with each other. We conclude that the differences 
between the three works are driven by systematic offsets 
in X-ray observables, and the slope of the f gas -M rela- 
tion. In particular, the statistical treatment of each of 
the samples — i.e. the fact that each data set used a dif- 
ferent statistical method for fitting the scaling relations 
and correcting for selection effects — appears to play a 
secondary role relative to the overall level of systematic 
offsets in the X-ray observables themselves. 

3.2. The M-Y x Relation 

The top panel in Figure [3] shows the 68% confidence 
region for the M-Yx relation of the three works con- 
sidered here, after subtracting a mean scaling relation 
defined by a self-similar slope of 3/5, and setting the 
amplitude to the mean amplitude of three works consid- 
ered here at Yx = 4 x 10 14 M keV. The parameters for 
this mean scaling relation are given in Tabic [TJ The V09 
and Pll-LS relations are in excellent agreement, despite 
the large differences in the cluster-by-cluster comparison 
of mass and Yx measurements, demonstrating that the 
systematic offsets from paper I move clusters along this 
best-fit relation. Unfortunately, because we do not know 

9 We note that when applying a shift in the mass-axis and a tilt, 
it matters whether wc shift first or tilt first. A uniform mass shift 
would only be appropriate for parallel relations, so we tilt first, 
then shift. 
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Fig. 3. — Top panel: 68% confidence band for the M— Yx relation 
derived from the Pll-LS, M10, and V09 data sets. For ease of 
comparison, we subtract out a reference mean relation as per the 
bottom panel in Figure [T] Bottom panel: As top panel, but after 
correcting for the systematic differences in cluster observables, and 
tilting the Mantz relation to that of an / gas <x M - 15 model. As 
in the case of the L-^—M relation, all scaling relations are now in 
good agreement. 



what is the source of the P11-LS-V09 offset in cluster ob- 
servables, it is difficult to guess why this "conspiracy of 
errors" would happen. As for the M10 relation, we once 
again see a significantly different slope, so the amplitude 
offset depends on the mass scale under consideration. 

The P11-LS-V09 comparison argues against instru- 
mental calibration being the principal source of the ob- 
served systematic errors. Specifically, a temperature 



bias, bx, induces a bia s 6m = b^ in hydrostatic masses 
(|Vikhlinin et al.l 120061 ). This in turn induces a bias 
6 oaB = b^f = b T 6 via aperture effects, so the net bias in 



Yx is by = b T . Letting a w 5/3 be the slope of the Yx~ 
M relation, we see that a temperature calibration bias bx 
will induce a bias b T ~ b T 9 in the amplitude of the 

Yx~M relation. In other words, an overall temperature 
bias due to instrumental calibration docs not appear to 
shift galaxy clusters along the scaling relation. 

The lower panel in Figure [3] illustrates the offset- 
corrected scaling relations. We have also tilted the M10 



Fig. 4. — Top panel: The Ygz-Yjc relation at z = 0.23 for 
each of the data sets we consider, as labelled. There are large 
differences between each of the data sets. Bottom panel: The 
Ygz— Yx relation after correcting for cluster observable systematic 
offsets, and correcting the Vikhlinin curve for the expected redshift 
evolution. The vertical dotted line is the pivot point of the V09 
relation. 

relation to that of a 7 = 0.15 / gas model. As in the 
case of Lx~M, after these corrections all three scaling 
relations are in agreement within the expected errors. 

3.3. The Ygz-Yx Scaling Relation 

The top panel in Figure |4] shows the Y^z-Yx scal- 
ing relation after subtracting out the reference relation 
(lnD 2 A Y sz \Yx) = ln(CFx) where C is given by equation 
and assuming both D A Ysz and CYx are measured in 
the same units. Because we assume a = 1 for the Mantz 
data, we only place a single point with error bars at the 
pivot point of the sample. There are large differences 
between the Pll-LS(z=0.23) and the V09 and M10 sam- 
ples. 

The bottom panel in Figure Q] corrects the V09 and 
MI0 scaling relations for the observable offsets from Ta- 
ble [U relative to the Pll-LS(z=0.23) data. The V09 scal- 
ing relation is modified using the shifts relative to the 
low-redshift PII-LS data, and then again by the red- 
shift correction between low and high redshift based on 
the Ysz/Yx ratio. The M10 value is shifted using the 
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Fx shift relative to Pll-LS(z=0.23) directly. With these 
corrections, all data sets agree on the amplitude of the 
Ysz/Yx relation, though some differences in the slope are 
apparent. 

3.4. Summary of Results 

The main result of this section is very clear: despite 
differences in fitting methods as well as differences in 
the treatment of selection effects, after accounting for 
the systematic differences in derived X-ray cluster ob- 
servables and modeling of the slope of the M gas -M rela- 
tion, all scaling relations arc in good agreement with each 
other. In other words, systematics in the treatment of 
selection effects and/or fitting methods are clearly sub- 
dominant to the overall systematic offset in the input 
data used to derive the scaling relations. 

4. THE Y sz -M RELATION AS CALIBRATED FROM 
X-RAYS 

Having specified the M—Yx and Ysz-Yx relation, the 
Ysz-M relation is nearly completely specified. We now 
derive expressions for the amplitude and slope of the 
Ysz~M scaling relation in terms of the others, which we 
then apply to each of the X-ray data sets considered in 
this work. 

Part of the motivation for this work is simply to trans- 
port the M-Yx calibration from X-ray data sets to Ysz - 
M, in order to enable cosmological interpretations of SZ 
cluster samples. Just as importantly, however, we wish to 
determine the differences in the predicted Ya^-M scaling 
relation for each of these data sets, as these differences 
will necessarily have an important impact in the inter- 
pretation of the SZ-richness scaling relations of maxBCG 
galaxy clusters (paper III). 



4.1. Method 



We define 



2/,= hi (Y x /Y ) 
y sz = ln(D 2 A Y sz /Al) 
m = In (M/M ) 



(6) 
(7) 

(8) 
(9) 

where the normalization constants Y , Aq, and M are 
chosen so as to decorrelate the amplitude and slope of 
the cluster scaling relations. In the previous section, we 
specified the Ysz-Yx and M-Yx scaling relation. Com- 
bining this two relations, we can arrive at an expression 
for the Ysz~M scaling relation. Specifically, let us write 

{y S z\m) = a sz \ m + a 8Z \ m m (10) 

{m\y x ) = a m \ x + a m \ x y x (11) 

(Vsz\yx) =a sz \ x + a sz \ x y x (12) 

We wish to solve for a sz \ m and a sz \ m in terms of known 
parameters. Using a local power-law model for cluster 
abun dan ces — se e Appen dix |A"] and specifi cally equa- 
tions |A5] and [All] (see also lWhite et al.l[2010| ) — one has 
that the parameters of the various relations above are 
related via 



^sz\x " 



^sz\m 



+ a 



'^sz\rn^m\x (-^) 



where r sz<x \ m is the correlation coefficient between y sz 
and y x at fixed to, j3 is the slope of the halo mass func- 
tion (dn/dhvM oc M - ^), and o m \ x and o~ m \ sz are the 
scatter in to at fixed y x and y sz respectively. Note the 
term in square brackets in equation [13] is simply the naive 
relation obtained from plugging equation QT] into equa- 
tion nm 

Since both a m \ x and a sz i x are known from the previous 
section, we can readily solve for a sz i m , 



^sz\x 
^m\x 



(15) 



To solve for a sz \ m from equation ll31 we must first be able 
to estimate the correction term depending on the correla- 
tion coefficient, which requi res that we know both r and 
°~m\sz- While equation IA13I in Appendix lAl relates these 
two quantities, one still needs an additional independent 
input to close the system. 

We arrive at an additional relation between r and o~ m \ sz 
using a simple model for h ow Ysz and Yx are c orrelated. 
Motivated by the fact that iRozo et all (2012aft find that 
the intrinsic scatter in the Ysz-Yx relation is consistent 
with zero, we make the assumption that the noise in y sz 
has two contributions: an intrinsic contribution Syi n t, 
which is identical and perfectly correlated with the in- 
trinsic noise in y x , reflecting variations in the details of 
the physical properties of the intra-cluster medium, and a 
an additional noise Syi ss that is due to structures along 
the line of sight to the clusters, and which is uncorre- 
cted with fluctuations in y x . Setting (Syf nt ) = c^,, 

(SyintSyus) = 0, Sy sz = Sy mt + Syi SSl and Sy x = 8y lnt , 
we find 

_ (Sy S zSy x ) _ o- x \ m , , 



sz.x \m 



x\m^ ' sz\m 



' sz\rn 



In addition, from equation IA6I in Appendix [Al we have 



' m\x 



(17) 



with a similar relation for y sz . Putting these two to- 
gether, we can rewrite equation [13] as 



^sz\m 



a 



s ~ i n LL rn x 



m\x ' 



(18) 



As a rough order of magnitude, we expect a slope of 
the halo mass function f3 w 3, and a m \ x « 3/5, which 
for 10% scatter in mass corresponds to a nearly negligi- 
ble 5% amplitude shift. Together, equations [TE\ and [15] 
completely specify the Y§z~M scaling relation in terms 
of known parameters (i.e. in terms of the Ysz-Yx and 
M-Yx scaling relation parameters). 

Let us now turn to t he sca tter in the Y$z~M scaling 
relation. From equation IA13I in Appendix [A] we have 



sz\x 



m\sz 



2r a 



s z , x | m ° in. | x ° m \ s z 



(19) 

Using equations [16] and IA6[ we can rewrite the r term 
above in terms of (J m \ x only. Solving for c^| sz , and using 



equation I A6I again to obtain o~ sz \ m from er r 
at 

2 2 2 ^sz\x {^sz\x 2) 

^ ~sz\m ® sz\x ® m\x ^2 



we arrive 



(20) 
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Fig. 5. — Likelihood distribution for the correlation coefficient 
r = Cov(y sz , yx)/& x \m a 'sz\m f° r each of the three data sets we 
consider, as labelled. The maximum likelihood value for r, as well 
as the upper and lower limits defining the 68% confidence contour, 
are summarized in Table [3] 



This gives cr sz \ m in terms of known quantit ies. To solve 
for r SZyX \ m , we again use equations \W\ and IA6I to write 



*sz\m 



in terms of : 



We arrive at 



and cr m \ x , and solve for r szx \ m . 

'1/2 



Q 



j | x 



sz\x ) 



(21) 



We evaluate the right hand side of the above equation 
using the central values from the V09 data set to arrive 
at T S z,x\m — 0.84. As expected, the fluctuations about 
the mean in Ysz and Yx at fixed mass are heavily corre- 
lated. The likelihood for the correlation coefficient r is 
shown in Figure [5j and the corresponding constraint is 
summarized in Table 

Having derived analytical relations between the input 
M-Yx and Ygz—Yx scaling relation parameters, and the 
Ysz~M scaling relation parameters, we now wish to prop- 
agate the uncertainty in these input scaling relations (see 
Table [1]) into uncertainty in the Ygz~M relation. The 
first step is to specify the choice of units for each data 
set. The mass units M are set to Af /10 14 M@ = 4.8, 
5.5, and 10.0 for the V09, Pll-LS, and M10 data sets. 
The slope of the halo mass function (3 is evaluated at 
this pivot point for our fiducial cosmology, and the cor- 
responding values are tabulated in Table [3] For the SZ 
units, we set A 2 , = 10 -5 Mpc 2 . We adopt the pivot point 
Yq quoted for the M-Yx relation. Our choice of units for 
Mq is set so as to match the pivot point in the Lx~M re- 
lation, unless doing so induces strong covariance between 
the amplitude and slope of the Y$z~M scaling relation. 
This only occurs for t he Pll-LS(z=0.23) da t a set, w here 
the mass pivo t of th e iPlanck Collaboration! (|2011b[ ) and 
lArnaud et al.l (|2010T l data are not well matched. In this 
case, we chose our mass units to decorrelate the ampli- 
tude and slope of the Y§z-M scaling relation, since this 
is the quantity of interest in this section. 

To determine the uncertainty in the Y$z~M scaling 
relation, we randomly generate the M-Yx and Ygz-Yx 
scaling relation parameters from the values quoted in Ta- 



ble [JJ These randomly drawn parameters are input into 
equations ([15]). (JTSJ and ([20) to solve for the Y sz -M 
scaling relation parameters. Table [JJ only quotes uncer- 
tainties for the scatter a m \ z for the M10 data set. For the 
iVikhlinin etHI (|2009al ) and I Arnaud "et~aTl (pOlOh M-Yx 
relation, we adopt a uniform prior on the variance cr 2 ,^ 

with <r^ x < 0.07 2 and a^ x < 0.09 2 respectively. As for 
the M10 Ysz~Yx scaling relation, we assume a conserva- 
tive flat prior on the variance c 2 . < 0.15 2 . We perform 

10 5 independent random draws to estimate the resulting 
likelihood distributions, also enforcing the prior \r\ < 1. 
The resulting Ygz-Af scaling relations are summarized in 
Tableland Figurc|Hl and the likelihood distributions arc 
shown in Figure[10] We emphasize that because the mass 
pivots arc different for different relations, one cannot use 
Figure [10] to determine whether the various Ysz~M rela- 
tion are consistent with one another. For that purpose, 
one should focus instead on Figure [H 

Finally, a warning. The relations between the ampli- 
tudes, slopes, and scatters of the various relations de- 
noted above are exactly correct. However, our Monte 
Carlo method for propagating the uncertainties in the 
Ysz-Yx and M-Yx scaling relations into the Y$z~M 
scaling relation explicitly assumes that the uncertainties 
of these input scaling relations are uncorrelated. How- 
ever, in all cases many of the same clusters were used as 
inputs in the analysis leading to the input scaling rela- 
tions were shared. So, for instance, some of the clusters 
that contribute to the Ysz~Yx calibration also contribute 
to the M-Yx calibration. In this case, the most correct 
analysis would be to do a multi-variate scaling relation 
analysis fit from the start, something which is not really 
feasible given the available data, and certainly beyond 
the scope of this work. With this caveat in mind, we 
proceed. 

4.2. Results 

Figurc[6]shows the 68% confidence contours of the Ysz - 
M scaling relations corresponding to each of the three 
sets of parameters summarized in Table El As usual, we 
subtract out a reference scaling relation, defined by a 
self-similar slope of 5/3. The amplitude of the reference 
scaling relation is defined by the mean amplitude at M = 
5 x 10 14 Mq between the V09, M10, and Pll-LS(z=0.23) 
data sets, and is summarized in Table[3] There is a clear 
amplitude offset between the V09 and Pll-LS(z=0.23) 
data which mirrors the difference in the Lx~M relation. 
The offset is partly caused by the systematic differences 
in mass calibration between Pll-LS and V09, and partly 
by the redshift evolution in Yx between low and high 
rcdshifts in the Pll-LS data. The M10 relation is in 
good agreement with the V09 relation where the samples 
overlap (M ~ 10 15 Mq), but the slope of the relation is 
very different, reflecting the differences in slope of the 
model. 

Comparing our self-consistently propagated scaling re- 
lation for the Pll-LS(z=0.23) data set to that published 
in Pll-LS, we see that ours has a higher amplitude, and a 
different slope. There are various reasons for this. First, 
an amplitude offset is expected given the difference in 
the Ysz/Yx ratio between low and high redshift cluster 
samples discussed in paper I. Likewise, a difference in the 
slope of the Ysz~Yx relation should also result in different 
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Fig. 6. — Top panel: 68% confidence regions for the Ygz— M scal- 
ing relation as determined by each of the four data sets we con- 
sider, as labeled. The p urple line is the Ysz-M scalin g relation 
at 2 = 0.23 as quoted in Planck Collaboration (2011b). Bottom 
pane l: Comparison of our derived scaling relations with the data 
from IMarrone et al.l 120111) . The purple lines show the 68% confi- 
dence regions of the M— Yjjz scaling relation from that work. Oth- 
erwise, the color key is the same as for the top panel: red=M10, 
blue=Pll-LS(z=0.23), black=V09. 



Ysz-M slopes. The role of selection effects is a bit more 
unclear: Pll-LS correct for selection effects by perform- 
ing Monte Carlo samples, which leads to a lowering of 
the amplitude by 7%. On the other hand, we estimated 
scatter corrections to the naive plug-in method — which 
we expect is roughly equivalent to directly fitting Ysz - 
M by setting all clusters to the mass M{Yx) — to lead 
to a 5% downwards correction, with some uncertainty 
associated with the value of the scatter. Because we ex- 
pect that the Monte Carlo procedure of Pll-LS would 
detect the scatter corrections introduced in this work, it 
is not surprising that the two corrections are nearly iden- 
tical, and it suggests that selection function effects that 
are over and beyond these scatter corrections are negligi- 
ble, as would be expected for scaling relations with tight 
scatter. 

In an attempt to shed further light on the Yqz~M scal- 
ing, the bottom panel of Figure [6] compares our predicted 
scaling relations to the LoCuSS sample measurements of 



iMarrone et al.l (|2011l) . That work e stimated M^qq usin g 
weak lensing shear measurements ([Okabe et al.l 120 j~0h . 
The individual data points are shown, while the solid 
purple lines span the 68% confidence region of the corre- 
sponding Ygz~M relation. Because the intrinsic scatter is 
small, the corrections from this inversion are completely 
negligible relative to the statistical errors. Consequently, 
we ignored the scatter corrections when estimating the 
Ysz-M relation parameters from the published M-Y^z 
relation parameters E3 It is clear from the figure that 
their best fit scaling relation is in agreement with all 
three of the scaling relations we computed, so the Lo- 
CuSS data is not yet of high enough quality to unam- 
biguously prefer one data set over another. 

In this context, we not e that a recent paper by 
iPlanck Collaborationl ()2012l ) has noted there are large 
systematic differences bet ween their X-ray m asses and 
the weak lensing masses of lOkabe et al.1 (|2010f ) . We treat 
this problem in more detail in paper III. Here, it suffices 
to note that these systematic differences are not sufficient 
to create statistic a l tens ion between our results and those 
of IMarrone et al.l (|2011l ) because of the large uncertain- 
ties in the latter work. 

5. THE Y sz -L x RELATION 

We now use the Y$z~M and Lx~M scaling relations 
to self-consistcntly derive the Y§z~Lx scaling relation. 
There are two motivations for this work. First, in the ab- 
sence of systematic errors, our prediction for the Ygz~Lx 
scaling relation must be fully self-consistent with the ob- 
served Ysz~£x relation. We have seen, however, that the 
scaling relations of galaxy clusters systematically vary 
from data set to data set. If we can determine that for a 
given data set our predicted Ysz^x relation is not con- 
sistent with the observed Y^z~Lx relation, then that is 
strong evidence that these exist systematic errors in the 
corresponding data set. 

To the extent that a multivariate Gaussian approxima- 
tion is a valid description of the scatter about mean scal- 
ings, any data set whose relations are not self-consistent 
(i.e., deviate from the expectations of Appendix |A"]) nec- 
essarily suffers from systematic uncertainties. By the 
same token, the Y"sz-£x scaling relation plays a crit- 
ical self-consistency constraint for the SZ-optical scal- 
ing relations of galaxy clusters, and will therefore pro- 
vide an important test of the interpretation of the 
IPlanck Collaborationl (|2011cD results presented in paper 
III. 

5.1. Method 
Similar to the previous section, we define 

l x =ln(L x /L Q ). (22) 

Having already specified the Y^z~M and the Lx~M 
scaling relations, the Y$z~Lx relation is constrained, up 
to the value of the covariancc between Ygz and Lx at 
fixed M. Because the scatter in Lx is dominated by the 

10 We did, however, implement a prior that the slope of the 
M—Ysz relation must be larger than a m | sz > 0.2. Because the 
amplitude of the Y^-M relation is — a m \szl a m\sz ( scc equation 
|A5}, allowing a m \ sz to get close to zero is both unphysical and 
numerically unsound. 
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TABLE 3 

Derived Cluster Scaling Relations at z ■ 



0.23 



Relation 



AO 



P 



In ip 



"'In,/'!) 



Sample 



D 2 A Y SZ -M 
D 2 A Y SZ -M 
D\Y SZ -M 
D 2 A Y SZ -M 



1.8 
5.5 
10.0 
5.0 



2.75 
2.93 
3.95 



1.34 ±0.07 
1.97 ±0.06 
2.54 ±0.20 
1.58 



1.61 ±0.11 
1.48 ±0.12 
1.48 ±0.09 

5/3 



0.12 ±0.03 
0.20 ±0.04 
0.15 ±0.03 



7^+ - 20 
u -'°-0.20 
f 0.15 
-0.20 
66+° 18 



0.60 H 



V09 

Pll-LS(z=0.23) 
M10 

Reference 



D 2 A Y SZ -L X 4.0 
D A Y SZ -L X 4.0 



DiY sz -L y 



10.0 



1.33 ±0.13 1.01 ±0.11 0.35 ±0.07 
1.45 ±0.15 0.92 ±0.10 0.32 ±0.10 
2.27 ±0.31 1.10 ±0.08 0.40 ± 0.08 



V09 

Pll-LS(z=0.23) 
M10 



a In all cases, we assume the V~X relation takes the form lni/> = In if>o + a ln(x/xo)- Our choice of 
units are 10 14 Mq for mass, 10 44 ergs/s for Lx, and 10 — 5 Mpc 2 for D a Yq z . All scaling relations 
arc appropriate for z = 0.23, the median redshift of the maxBCG cluster sample. The quantity f} 
is the slope of the halo mass function at the pivot scale of the Y% z —M relation. All errors are the 
standard deviation of the distribution, except for the correlation coefficient r between Y"sz and Vx 
at fixed mass. The likelihood distribution for r is highly non-gaussian (see Figure 0, so we quote 
the maximum likelihood value, and the error bars define the 68% likelihood contour. 



details of the X-ray emission from the core of galaxy clus- 
ters, while the scatter in Ysz is dominated by projection 
effects, the two are not likely to be correlated at any sig- 



nificant level, so one might expect r s 



,l\r 



0. On the 



other hand, iStanek et al.1 (|2010T l derived theoretical pre- 
dictions for the correlation coefficient between Ysz and 
Lf,oi, the bolometric X-ray luminosity, based on numer- 
ical simulations. They find r « 0.8, which reflect the 
fact that clusters with more gas and/or hotter gas will 
have both higher X-ray luminosity and higher pressure. 
In the absence of more detailed information, we adopt a 
uniform prior r sz j\ m G [0, 1]. 

Using the for mulae in Appendix |A] (sec also 
IWhite et al.l l2010h . and following arguments similar to 
those in the previous section, we can write the ampli- 
tude, slope, and scatter, of the Ysz~Lx scaling relation 
in terms of the parameters characterizing the Y^-M and 
the Lx~M scaling relation. In the notation of the previ- 
ous section, we find 



a sz\l : 



Usz\l ■ 



C*sz\m 



, 2 



a sz\l — a sz\f 



<? m \i [r SZt i\ (23) 

(24) 

(25) 



a m\l + a m\sz ~ 2 r szd\mVm\lVm\sz 



The term in square brackets is the naive relation one 
would obtain in the absence of scatter, and the scatter 
in mass at fixed observable is relat ed t o the scatter in 
observable at fixed mass via equation lA6l To get an order 
of magnitude estimate for the correction term, we set 
= 3, a. sz \ m = 5/3, r szAm = 0, and a m \i = 0.25, which 
yields a correction term ps —0.3, corresponding to a 30% 
offset in Ysz at fixed Lx, a very significant correction. If 



we set r. 



,l\r 



1 instead, and adopt a n 



0.1, the 



correction is only reduced to « —0.2, still large. 

To apply this formalism, we must again specify units. 
We adopt as mass units the pivot point of the Ysz~M re- 
lation, which matches that of the Lx~M relation for the 
V09 and M10 data sets. For the Pll-LS data set, the 
pivot point of the Lx~M relation is lower than that of 



the Ysz~M relation. When sampling the Lx~M scaling 
relation, we randomly draw the amplitude and slope at 
the pivot point of the Lx~M relation, and then shift the 
amplitude to the new pivot point. This procedure explic- 
itly includes the induced covariance between amplitude 
and slope at the new pivot point. The scaling relations 
are sampled 10 times, and the corresponding Ysz~Lx 
parameters constrained from the Monte Carlo distribu- 
tion. The units for Lx are chosen so as to decorrelate the 
amplitude and slope of the Ygz~Lx relation, and are set 
to L Q = 3 x 10 44 ergs/s for the V09 and Pll-LS data sets, 
and Lx = 10 45 ergs/s for the M10 data set. The units for 
D\Ysz are set to Aq = 10 -5 Mpc 2 in all cases. As in the 
case of the Ygz~M relation, we note that in principle the 
uncertainties in the input scaling relation parameters can 
be correlated, for instance due to shared clusters in the 
calibration of the Lx~M and Fgz-M scaling relations, so 
the same caveat mentioned in the previous section holds 
for this discussion as well. 

5.2. Results 

Figure [TT] shows the likelihood for each of the data 
sets we consider. There are significant differences in 
amplitude between the various works, but the uncer- 
tainties are al so large. We compare each of our pre- 
dictions to the iPlanck Collaboration! (1201 lbfl data. Fig- 
ure [7] s hows all galaxy clusters in IPlanck Collaboration! 
(|2011bD in the redshift slice [0.13,0.3]. We adopt this 
cut to ensure that we only look at clusters near the red- 
shift of interest (z = 0.23), and also because from paper 
I w e know that the low and high redshift populations of 
the lPlanck Collaboration! (|2011bf ) clusters appear to have 
different properties. Note that because Ysz depends on 
i?50o, when comparing our predictions to the Pll-LS Ysz 
data we need to account for the systematic difference in 
mass between the various groups as per the Appendix B 
in paper I. We do so by shifting the data points rather 
than our predictions, though we also note that because 
Ysz depends only mildly on i?5oo, the shifts due to aper- 
ture corrections are much smaller than the differences 
between the predicted scaling relations (i.e. aperture cor- 
rections in Ysz are a second order effect). 

In addition to showing the data points themselves, Fig- 
ure [Jj shows as a dark-green line our best-fit relation to 
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Fig. 7. — Comparison between the Planck Collaboration (2011b) 
data for clusters in the redshift range z £ [0.12,0.3] and the Y"sz- 
Lx relations derived sclf-consistcntly from t he Y"sz-M and Lyj-M 
relation s. The green line is our best fit to the Planck Collaboration 
Il2011bl ) data, while the black solid (dotted) lines mark the la (2a) 
confidence regions of the mean scaling relation. Blue points are 
cool-core systems, while red points are not. 

the Pll-LS(z=0.23) data. The amplitude of the fit is cor- 
rected for selection effects based on the results of Pll-LS, 
which amount to a 9% reduction of the best fit amplitude. 
The la and 2a regions for each of the three predicted 
scaling relations are shown in Figure [7] as black solid and 
dotted lines respectively. The best fit relation to the data 
is clearly that derived from the Pll-LS(z=0.23) scaling 
relations, but the best fit relation is within the 2a regions 
of the predicted scaling relations for all three data sets. 



Consequently, we are unable to identify which data sets 
suffer from systematic errors using the Ys z~Lx measure- 
ments from lPlanck CollaborationT^OllbO . 

We have also attempted to compare our self- 
consistently propagated scalin g relations to the data from 
IPlanck Collaboration! (|2011aD . henceforth referred to as 
Pll-X. Pll-X assembled a heterogenous X-ray c luster 
catalog, the MCXC catalog (jPiffaretti et al.ll201lD . and 
used it to measure the Y§z~Lx relation. Unfortunately, 
this comparison is not trivial. As noted in that work, 
the Ysz~Lx scaling relation from the MCXC sample and 
the Pll-LS scaling relations are not consistent, with the 
MCXC relation shifted towards more luminous systems 
(lower amplitude). This is illustrated in Figure which 
show the 2a regions of the Pll-X results as the solid pur- 
ple lines O Also shown as purple squares with errors are 
the binned data from the same work. The dashed dark- 
green lines are the 2a confidence regions as reported in 
Table 2 of IPlanck CollaboratiorJ (|2011bl . i.e. the Pll- 
LS result). The z 6 [0.12,0.3] clusters in Pll-LS arc 
shown as circles with error bars, and our own fit to the 
data is similar to that of Pll-LS. We see that the 2a 
bound aries of the two relations — t hat of Pll-LS and 
that of IPlanck Collaborationl (|2011a| ) — barely overlap. 
The amplitude difference at the pivot scale of the Pll- 
LS data (L x = 7 x 10 44 M s ) is AlnYgz = 0.29, a 3.2er 
offset. 

Pll-X attributes the offset in their result relative to 
that of Pll-LS to selection effects in the MCXC sam- 
ple. However, we believe this interpretation is not cor- 
rect for several reasons. First, the Pll-X sample relied 
on archival data, which heavily favors X-ray selected sys- 
tems. Thus, it is unclear whether the MCXC catalog is 
really any more "X-ray selected" than the Pll-LS sam- 
ple. More importantly, we emphasize that we are con- 
sidering the Y$z~Lx relation: that is, Ysz is consider the 
dependent variable. By definition, if our sample contains 
all galaxy clusters of a given X-ray luminosity, we can ex- 
plicitly compute {Yq,z\Lx) without any corrections, i.e. 
there are no selection effects of which to speak of. Finally, 
Pll-X argue that the reason for the offset is a larger pres- 
ence of cool-core systems in their data relative to Pll-LS. 
We have refit the z € [0.13, 0.3] data using only cool-core 
systems. The re sulting best fit relation is sti ll slightly 
to the left of the IPlanck Coi laborationl (|2011a[) relation, 
with an amplitude offset Am Y" S z = 0.08 ± 0.14. This is 
consistent with zero, but t aken at face value, wou l d argue 
that all galaxy clusters in IPlanck Collaborationl (|2011a| ) 
must be cool-core systems in order to account for the 
observed offset. 

Regardless of these difficulties, we can still compare the 
Pll-X scaling relation to our predicted V09, M10, and 
Pll-LS(z=0.23) Ysz-Lx relations. As a simple illustra- 
tive example, Figure [8] shows as a black line the V09 
prediction for the Y$z~Lx scaling relation. In all cases, 
we find that the Pll-X result is in better agreement with 
our predictions than the raw Pll-LS data. 

In short, it is not clear what drives the difference be- 
tween the Ysz~Lx relations of Pll-X and Pll-LS, nor 

11 To compute these regions, we used the fits from Table 4 in 
Pll-X where both the amplitude and slope of the relation are al- 
lowed to vary. We also corrected to account for the difference 
between (lsz|-^x) — the quantity measured by Pll-X — and 
(lnYgzlix). 
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Fig. 8. — Comparison between the Pll-LS (dark-green dashed 
lines) and Pll-X (purple solid lines) fits for the Ygz—Lx relation 
at z = 0.23. Both bands represent 2<r confidence intervals. Blue 
points are cool-core systems, while red points are not. At 2<r, the 
two regions barely overlap. The difference in amplitude at the pivot 
point of the Pll-LS data (L x = 7x 10 44 ergs/s) is AlnVsz = 0.29, 
which is significant at the 3.2cr level. The black solid line is the 
predicted Ysz-Lx relation from the V09 scaling relations, and is 
indicative of the excellent agreement between the predicted Ysz- 
Lx relations and the Pll-X results. The agreement of Pll-X with 
M10 is also excellent. 

which Ygz-Lx scaling relations arc ultimately correct. 
The best we can say is that all three predicted Ysz _ Lx 
scaling relations — V09, M10, and Pll-LS(z=0.23) - 
are viable, with better agreement found between our pre- 
dictions and the Pll-X results. We will return to the 
problem of the Ya^-Lx relation in paper III. 

6. COSMOLOGICAL CONSISTENCY 

Having performed a detailed comparison of the clus- 
ter scaling relations from three different data sets, we 
test which of these relations are also consistent with the 
currently favored WMAP7 ACDM cosmological model. 
This t est is not a priori : becau se of the cosmological ana l- 
ysis in lVikhlinin et al.l (|2009b[ ) and lMantz et alJ ()2010aD . 
we already know that both the V09 and M10 scaling re- 
lations are consistent with WMAP7. Here, we confirm 
this agree ment, and t e st the additional Lx~M scaling 
relation of iPratt et "ail (|2009D . which uses the same mass 
calibration as Pll-LS. 

We proceed as follows. We randomly sample the 
CMB +BAO+ffn Monte Ca rlo Markov Chains (MCMC) 
from iKomatsu et~aTl (|2011l ) for a "vanilla" flat ACDM 
cosmology. The parameters that we vary are as, n s , Q m , 
h, and fib- For eac h random sampling, we compute the 
iTinker et all ([20081 ) mass function. We also randomly 
sample the amplitude, slope, and scatter of the Lx~M 
relation for each of the V09, M10, and Pll-LS scaling 
relations, and convolve the resulting mass function and 
P{Lx\M) relations to arrive at the predicted X-ray lu- 
minosity function, 

The uncertainty in the predicted X-ray luminosity func- 
tion is estimated from the variance over our Monte Carlo 
realizations. 



We compare this prediction to the X-ray lumi- 
nosity function from the REFLEX cluster survey 
(|Bohringer et al.l 120021) . Note that both the reported 
P(Lx\M) scaling relations and the REFLEX luminos- 
ity function measurement assume f2 m = 0.3 rather than 
the WMAP7 central value tt m = 0.27. We do not ex- 
pect this difference to be an important systematic. For 
instance, the difference in the luminosity distance out to 
z = 0.23 between the two models is only sa 0.5%, so we 
do not expect more than a few percent changes over the 
range of cosmologies considered here. These differences 
are much smaller than the relative differences between 
the Lx~M relations considered here. 

Figure [9] shows the comparison between our predicted 
X-ray luminosity functions, and the REFLEX luminosity 
function. For our predictions, we evaluate both the mass 
function and Lx~M scaling relations at z = 0.08, the 
median redshift of the REFLEX clusters. The depen- 
dence of our results on the assumed redshift is mild. The 
bands for our theoretical predictions show the 68% con- 
fidence intervals, estimated from 10 3 random samplings 
of the cosmological and scaling relation parameters. The 
width of the REFLEX band is computed by random sam- 
pling of the Schechter luminosity functio n parameters 
and a using the quoted uncertainties in iBohringc r et al.l 
(|2002D . In addition, it is important to note that while 
in our work Lx was defined using an aperture -R50O) the 
REFLEX catalog defines Lx using growth curve analy- 
sis, so the two luminosities are not directly comparable. 
Because Lx varies very slowly with aperture, we expect 
that a ±10% systematic uncertainty in the luminosity is 
a reasonable estimate the impact that this difference can 
have. We convert this uncertainty into an uncertainty 
in the luminosity function using the logarithmic deriva- 
tive of the average luminosity function, and add this in 
quadrature to the statistical error to arrive a the total 
uncertainty in the empirical luminosity function. 

The small vertical lines along the cc-axis denote the 
Lx corresponding to the mass pivot point of each of the 
three Lx~M relations. For properly calibrated scaling re- 
lations, the predicted abundance function should be sta- 
tistically consistent the observed abundance at the pivot 
point. As expected, the V09 and M10 models are in good 
agreement with WMAP7, with the observed offsets sig- 
nificant at 0.6cr and 0.3c respectively. The Pll-LS model 
is in mild tension (2.1ct) with WMAP7. 

We have repeated this experiment using the addi- 
tional cosmological priors from the Bary on Oscillation 
Spect roscopic Survey (BOSS) experiment (| Sanchez et al.1 
l2012f) . Because the full likelihood is not yet published, 

1 /2 

we sample only the cosmological parameter a^flm , cor- 
responding to the standard cluster normalization con- 
dition. The cosmological prior from BOSS for this pa- 
1/2 

rameter is cF^m = 0.441 ± 0.013 (A. Sanchez, private 
communication). All other cosmological parameters are 
held fixed to their fiducial value. We have further veri- 
fied that when we perform this simplified analysis on the 
WMAP7 results, our results are very nearly identical to 
the results when we allow all cosmological parameters 
to vary. The addition of BOSS data sharpens the pre- 
vious discussion: V09 and M10 remain consistent with 
WMAP7+BOSS, with the abundance offset being signif- 
icant at 1.5cr and 0.7(7 respectively. The offset relative 



14 



Rozo et al. 



10" 6 



* V 1 1 1 1 


1 ' 




vuy 




Pll-LS 




M10 ~ 




REFLEX : 


r 


















\\ vV\ "v. 






z=0.08 










\ \ \ x 










. . j 1 


i ^ \ N 


1 


10 



L x (10 44 ergs/s) 

Fig. 9. — Comparison of the X-ray luminosity predicted from 
the Lx~ M relation derived from each of the data sets to the RE- 
FLEX luminosity function, as labelled. The bands shown corre- 
spond t o 68% confidence, and predictions assume a WMAP7 cos- 
mology (Komatsu ct al. 2011) sampled over erg, f2 m , h, n a , and f2i,. 
The vertical lines along the x-axis mark the pivot point of each Lx_- 
M relation. Both V09 and M10 are consistent with cosmological 
expectations, while the Pll-LS relation is in mild tension (2.1ct) 
with WMAP7. Adding BOSS priors (see text) increases the tension 
with Pll-LS to 3.3<t, while both V09 and M10 remain consistent 
with the REFLEX luminosity function. Both the mass function 
and the scaling relations are evaluated at z = 0.08, the median 
redshift of the REFLEX catalog. 



to the Pll-LS prediction is 3.3er, in modest tension with 
the WMAP7+BOSS results. 

Interestingly, an examination of Figurc[T](Lx _ Af), Fig- 
ure ED (M-y x ), Figure H (y sz -F x ), Figure E (Y sz -M), 
and Figure [7] (Fsz-ix), reveal that at the pivot point of 
the M10 data set, the amplitude of the V09 scaling re- 
lations are in good agreement with the amplitude of the 
M10 scaling relations. That is, the V09 and M10 rela- 
tions are in good agreement for the most massive clusters, 
but diverge at low masses because of the differences in 
the fgns-M relation between the two works. From Figure 
[H we see that the f gas model of V09 is a better match to 
the full set of scales probe d by REFLEX. Th is stands is 
in contrast to the results of lAllen et al.l (|2008l ). However, 
the latter work focused on the most massive galaxy clus- 
ters in the Universe, so the two results may be reconciled 
if there is a flattening of the slope in the f gas -M relation 
at the highest masses. 

7. SUMMARY AND CONCLUSIONS 

We have performed a detailed comparison of the 
X-ray scaling relations f r om t hree different data 
sets: Pll-LS dPratt et al.1 120091 lArnaud et all l2010l: 
Planck Collaboration! I2011bl)~ V 09 (iVikhlinin et al.1 
2009aHRozo et al.ll2012al ). and M10 (jMantz et al.ll2010bO . 
We compare the Lx~M, M-Yx, and Fsz~^x scaling re- 
lations and find varying degrees of tension between the 
works. Differences are partly traced to the systematic 
offsets in cluster observables characterized in paper I, 
but the constant / gas model of M10 drives slope differ- 
ences relative to V09 and Pll-LS. After correcting for 
these two effects, all cluster scaling relations are in good 
agreement with each other. Indeed, at the pivot point of 
the M10 scaling relations, all of the V09 and M10 scal- 



ing relations are in good agreement with each other: the 
difference between these two works is pronounced only if 
one extrapolates the M10 scaling relations to low mass 
systems. 

Having identified the sources of tension, we use these 
scaling relations to self-consistently recover the Y§z~M 
and Ygz—Lx scaling relations within the context of a 
multivariate Gauss ian model for clus ter properties (see 
Appendix El also IWhite et al.l |2010T) . In the limit of 
very small scatter, the transfer relations are equivalent 
to "plugging-in" one relation into another. In the pres- 
ence of finite scatter, there are additional terms involving 
property covariance and the local slope of the mass func- 
tion that must be taken into account. For instance, in 
the case of the Y$z~Lx scaling relation, these corrections 
arc of order 25%. 

The differences in mass calibration between the various 
data sets, presented in paper I, lead to different predic- 
tions for the Fgz-M scaling rela tion. Independent res ults 
from the LoCuSS collaboration ([Marrone are 
not yet able to distinguish between the various predic- 
tions. 

For Ysz-Lx scaling, we find that the differences be- 
tween the various data sets are moderate. Indeed, we 
compared each of the predicted Ygz~Lx scaling relations 
to the Pll-LS data, and found reasonable agreement in 
all cases. We also compar ed our predicted Yrz-Ly scal- 
ing relations with those of iPlanck Collaborationl (j2011al 
referred to as Pll-X), noting that there are large differ- 
ences between the Pll-LS and the Pll-X results ^.In- 
significance), finding much better agreement with the 
Pll-X results for the V09 and M10 data sets. The selec- 
tion effects required to reconcile the Pll-LS and Pll-X 
would have to be large, with the MCXC catalog com- 
prised exclusively of cool-core clusters, or SZ selection 
leading to an observed 30% amplitude shift; either solu- 
tion seems too large to be plausible, as correction terms 
from selection effects in SZ typically scale as the variance 
cr^ z | TO w 0.01, where as the observed offset is nearly 30% 

in Y S z- 

Finally, we consider whether the scaling relations from 
each of the three data sets are consistent with cosmo- 
logical expectations. We convolve the three Lx~M re- 
lations with the mass function predicted by a WMAP7 
cosmology — sampling over cts, fi m , n s , h, and Qb ~ 
and compare these predictio ns against the observed RE- 
FLEX luminosity function (Bohrin ger et al.ll2002f) . As 
expected, the M10 and V09 scaling relations are consis- 
tent with the REFLEX luminosity function at the pivot 
point of the two samples. The prediction from the Pll- 
LS Lx~M scaling relation is offset from the REFLEX 
data by 2.1er. Adding BOSS priors increases the tension 
with Pll-LS to 3.3cr, while both the V09 and M10 pre- 
dictions remain consistent with the data. In addition, 
comparing the M10 and V09 predictions, it is clear that 
the abundance of low mass systems is correctly predicted 
by V09 but not by M10, which argues that an / gas model 
/g as oc M 7 with 7 0.10 — 0.15 is a better match to the 
data over the full range of scales probed by REFLEX 
galaxy clusters. Because f^* is o bserved to be constant 
at high masses (jAllen et al.ll2008l ). this may be signaling 
a steepening of the f gas ~M relation at low masses. 

The fact that the three different X-ray data sets con- 
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sidered here and in paper I exhibit systematic differ- 
ences in mass calibration and cluster scaling relations 
has important consequences for the predictions of the 
YS7-N2CH) scaling rela t ion for optical galaxy clusters in 
iPlanck Collaboration! (|2011d ). In particular, the fact 
that the different data sets give rise to different Ygz~ 
M scaling relations also implies that the predicted Y"sz~ 
./V200 relations from each of these data sets will be differ- 
ent. In paper III, we explore whether any of these three 
data sets can resolve this issue, and whether doing so re- 
sults in a self-consistent picture of multi-variate cluster 
scaling relations. More specifically, we will will demand 
not just that the predicted Ysz-A^oo relation match ob- 
servations, but also that any two scaling relations can 
be combined to successfully predict the third, while si- 
multaneously satisfying cosmological expectations for the 
counts of galaxy clusters selected by any property. 
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APPENDIX 

A LOCAL MODEL FOR MULTI-VARIATE SCALING RELATIONS 

We consider the problem of a locally power-law halo mass function dn/d\aM oc M _/3 , and a vector of observable 
signals S, e.g . S = \Lx - Yx, Mp- a . K , Y r z, etc. }. The results summarized here are an extension of those presented in 
IStanek et al.l (|2010D and lAllen et al.l (| 2 1 If) , and are basically equivalent to those in Appendix C of iWhite et al.l 
(|2010T l . We define [i = ln(M/M ) and s = \n(Si/ Si, re f) where M and Si, re f represent a choice of units. The scaling 
relations between S and M are governed by the probability distribution P(s\m) which we assume to be Gaussian. The 
means of the distributions are parameterized as 

{s\n) = a + a /.i. (Al) 

The scatter is characterized by a covariance matrix, C, which has the property variance, af, along the diagonals and 
off-diagonal terms, Cy = ((sj — (si))(sj — (sj))). The correlation coefficient, given by = Cy/<XjCTj, lies between —1 
and 1. 
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Using Baycs Theorem, we can relate the S-M scaling relation to the MS scaling relation 

_ P( B \n)P(ti) P(s| M )P(/x) , 

{m ~ p(s) - j> pWW)- { } 

For a locally power-law model, P(/j) oc exp(— /3/x), the resulting probability density is Gaussian with mean and variance 



( M |s) = aC- 1 (s-aK |s -^ 



^(aC- 1 *)" 
For a single property, s, the above expressions reduce to 

{m\s) = 



(A3) 
(A4) 

(A5) 
(A6) 



The correction term to the mean mass, —fio , is the standard Malmquist bias correction. Here, is the scatter in 
halo mass of a sample selected by property s. 
The space density of clusters with multiple properties, s, is 



n(s) = dfi P(s\(i)n((i) 



Act,, 



( 27r )(JV-i)/2| C |i/2 
In the case of a single signal s, this reduces to 



exp 



(s-a)C-\s-a) 



n(s) = — exp 
a 



s — a 1 
2 



(A7) 



(A8) 



Note the slope of the abundance function is simply /3/a, as we would expect. There is a correction term to the 
amplitude, /3 2 c 2 , which reflects the number boost due to low mass halos scattered upward. 

Of particular interest to us is the case in which there are two signals Sx and s 2 . Let a^x and cr M 2 be the mass scatter 
in each, respectively. Then the mean selected mass, (/x|s) , and the mass variance for joint signals, s\ and s 2 , arc 



ajotx 1 (s 1 - ai) + a Ja 2 1 (s 2 - a 2 ) - ra Jo J(a t 1 (s 1 -ax) + a 2 1 (s 2 - a 2 )) 
(pt\s\,s 2 ) = 



1-r 2 



Finally, a further application of Bayes theorem allows us estimate P(s 2 \sx) 

1 



P(s 2 \sx) = P(sx, S2 )-J— = f dm P( Sl ,s 2 \fi)P(pL) 

P[Sx) P(Sx) J 



P(sx) P(sx) 

from which we find that P(s 2 \sx) is a Gaussian distribution with mean and variance 

(s 2 \sx) = a 2 \ m + a 2 ((/i|si) + r^o^^a) 
°2|i = a l + °\ 2 - Zro^xv^] 
The correlation coefficient between s 2 and m at fixed sx is given by 



r 2,u|l 



[l-r 2 + KiK2-r) 2 ] 1/2 ' 



(A9) 
(A10) 

(All) 

(A12) 
(A13) 

(A14) 



LIKELIHOOD PLOTS 

Figures [TBI and HT1 show the posterior distribution for the Ygz~M and Ygz~Lx scaling relation parameters summarized 
in Table O 
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Fig. 10. — Likelihood distribution and likelihood contours for the parameters of the Ygg— M scaling relation. Black solid=V09, blue 
dash-dot=Pll-LS(z=0.23), and red dashcd=M10. All contours arc 68% and 95% confidence. Note the different scaling relation assume 
different pivot points, so one cannot directly compare the amplitudes quoted between all data sets. 
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Fig. 11. — As Figure ITOl but for the Ygz— Lx parameters. 



